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I. TORQUE-SPEED CURVE MEASUREMENT 

o 

The measurement of the torque-speed curve is usually done by fixing the cell to a glass 
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slide and tethering a polystyrene bead to the flagellar hook. An optical trap monitor the 
rotational speed of the bead and the motor torque is calculated from r = (£ L +C,r)u ~ Cl^, 
where £r is the drag coefficient due to the internal friction in the motor, £l is the bead drag 
coefficient, and uj is the angular velocity. The torque-speed curve is then obtained changing 
£l by varying the bed size |l| or the viscosity of the external medium. An alternative method 
is to tether a cell to a glass coverslip by one of its shortened flagellar filament and expose 
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the cell to a rotating electric field [2]. Then the motor is broken spinning the cell backward. 
The difference between the cell body speeds before and after the broke of the motor, at the 



same value for the external applied torque, is proportional to the motor torque. 
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II. HOOK SPRING COMPLIANCE 



o . 

Q\ , Fig. SI shows the compliance curve F (motor torque versus angular displacement 9 — 9l) 



similar to the experimental one [3]. The following expression is used in our simulations: 



F(x) = kqx + Ax(ki — K )log[l + exp{(x — x )/Ax^], (1) 

with x = 6 — 9 L , k the spring constant a low load, /«i(3> «o) the spring constant a high 
load. At both low and high loads the spring behaves linearly with different spring constants 
K and K\ respectively. Ax is the angular displacement interval of the non linear region 
centered around the turning point xq. For the values of these parameters see Tab. SI. 

III. DISTRIBUTION FUNCTIONS OF t m AND t w 

Fig. S2 shows the distribution functions for t m and t w at different values of the load for the 
case N=l. Fig. S2(a) shows the average waiting and moving times. The arrows labeled with 
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the letter (b), (c), and (d) indicate the points (i.e. the speed values) where the probability 
distributions for t m and t w are shown in Fig. S2 (b), (c), and (d) respectively. The averaged 
waiting-time decreases slightly with the speed, while the averaged moving-time decreases by 
four orders of magnitude. The two time scales crossover at a speed around 170Hz, which 
naturally defines two regimes: i) (t w ) <C (t m ) and ii) (t w ) ^> (t m ). The crossover speed 
corresponds roughly to u n . 

Figs. S2(b), (c), and (d) show the differences in the distribution functions for t rn and t w 
at high, medium, and low load respectively. The waiting times are exponentially distributed 
because the waiting time interval is determined by independent chemical transitions, i.e., 
by Poisson processes. The average waiting time thus depends on the stator jump rate k, 
which varies between two constants k + and /c_ (except for extreme high load where k — 0): 
(t w ) oc 1/k. This explains why the averaged waiting time only weakly depends on the load. 
At low load, premature jumps are rare and the average angular movement S m is determined 
by a single stator jump and has a peaked distribution centered around Sq/N. This explains 
the peaked distribution for t m at low load, as shown in Fig. S2(b). At medium load, both 
(t m ) and the value of t m corresponding to the peak of the distribution increase, as shown in 
Fig. S2(c). At high load, the t m distribution develops a flat region for shorter time intervals, 
as shown in Fig. S2(b). This is a consequence of the decreasing slope of the total potential 
felt by the rotor on the positive force side. As a result, fluctuations of the rotor angle 9 due 
to thermal noise increase, which leads to many short moving time intervals. 

IV. DEPENDENCE OF THE TORQUE PLATEAU REGION ON K+ AND 5 C 

In order to understand the origin of the torque plateau region, we have studied the 
dependence of uo n on the ratio r = k + /k_ and on the cutoff S c for the case N=8 (similar 
results have been obtained for different values of N). We define u n as the speed value at 
which the torque decreases 10% from its value at stall. The size of the plateau regime is 
characterized by the following quantity: 

£ = ^n/^max- (2) 

Fig. S3 (a) shows the torque-speed curves for N= 8 and for two different values of r: r = 
and r = 1.2. For r = the torque, after a small plateau, decreases linearly with the speed. 



2 



Increasing r, u n increases and the torque-speed curve increase its concavity. In Fig. S3(b), 
E increases from 0.3 to 0.6 as r increases from 0.05 to 1. This corresponds to the values of r 
that can produce a well defined shoulder and at the same time maintain the independence 
of u> max on N, i.e. |A| < 1 (see Fig. 4(b) in the main text). 

Fig. S3(c) shows four torque-speed curves for different values of the cutoff S c . Starting 
from the value used in the main text, i.e. S c = So, uj n increases with S c . In particular, E 
increases significantly from S c = So to S c = 3So, at which the plateau size (ou n ) reaches a 
value slightly bigger than 60% of the maximum speed (Fig. S3(d)). 

In conclusion the extension of the torque plateau region increases with the value of rate 
k + and the cutoff S c , in consistent with our theory. 



V. ROBUSTNESS OF THE RESULTS AGAINST DIFFERENT ROTOR-STATOR 
POTENTIALS AND LOAD-ROTOR FORCES 

In order to verify the independence of our results on the specifics of the rotor-stator 
potential, we studied our model with asymmetric potentials and a smoothed symmetric 
potential with a parabolic bottom (instead of the V-shaped bottom). In the asymmetric 
case, the slope r+ of the left branch of the potential is much smaller than the slope r_ of 



the right branch, similar to the potential used in [4] (see the values used for r + and r_ in 
Tab. SI). Our model with the asymmetric potential yields qualitatively similar results as 
with the symmetric potential used in the main text. In particular, the maximum speeds 
near zero load are independent of the number of stators (see Fig. S4) provided the stator 
jumping rates satisfy: ^> t_/t + . Such a requirement can be understood intuitively 

in the following way. Given the condition r_ ^> r+, the force equilibrium (waiting phase) is 
achieved by having one stator spending part of its time dragging the rotor while all the other 
stators are pulling the rotor. The waiting period ends when this dragging stator jumps with 
a rate that depends on the fractions of time it spends on the two sides of the potential, which 
depend on the ratio t_/t + . Therefore, the condition that the maximum speed is dominated 
by k- (instead of k + ) has to be weighted by the ratio t_/t + . 

We have also studied a "semi-parabolic" potential V P (A9 = Q — Ql) (see insert in Fig. S5) 
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defined as: 

vP \ r [\A6\ - (Aflo/4)] if \A9\ > A6 /2, 
) t A6 2 /A6 if \A6\ < A6 /2, 

where tq is the positive slope of the symmetric potential, and A$o is the angular interval 
of the parabolic region centered around the bottom of the potential. Correspondingly, the 
chemical rate is a continuum function of AO: 

k{A6) = k + + (k- - k+)/[l + exp( - (A0)/A0 o )] (4) 

Fig. S5 shows the torque-speed curves for N = 1, 2, 8. The curves show the same charac- 
teristics as for the V-shaped potential shown in the main text. The plateau region is a little 
wider. This is due to the change of the slope near the potential bottom. A lower value of 
the slope slows down the motion of rotor, increasing the premature jump probability before 
it reaches the bottom. 

Next, we considered different forms of the force function F between the load and the 
rotor. Fig. S6 shows torque-speed curves for two cases: with and without a spring between 
the load and the rotor for different ratio k+/k~. The case with spring between the load and 
the rotor is studied in the main text; the case without spring corresponds to infinite spring 
constant (rigid connection between load and rotor) with the following equation for the rotor: 

d6 1 dV 



dt U + £ L d6 



+ \/2k B T/(S R + £ L )a(t). (5) 



Contrary to the model proposed in J4], the concavity of the torque-speed curve does not 
depends on the strength of the hook spring: the torque-curves are almost identical with and 
without spring. Instead, the concavity of the torque-speed curve depends on the ratio of 
the jump rates k + /k_, and also on the cutoff 5 C as shown before in Fig. S3. In particular, 
as shown in Fig. S6, for k + = the concavity is zero, and it increases as the ratio k + /k^ 
increases. 
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Torque (pN-nm) 



FIG. 1: Compliance curve. The graph shows the angular displacement, — 6l , between the rotor 
and the load as a function of the rotor torque. The non linear spring behavior follows approximately 
the experimental measurement in reference The curve corresponds to the function F(x) — F(0) 
used in our simulations. 
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FIG. 2: Waiting-time and moving-time statistic, (a) The average waiting-time (dashed line) and 
moving-time (solid line) averaged over 500 revolutions as a function of the rotational speed for 
N= 1. With increasing speed, the averaged waiting-time slowly decreases from 1ms to 0.1ms while 
the average moving-time decreases much faster over four orders of magnitude, from roughly 50ms 
to 0.005ms. The vertical arrows labeled with the letters (b), (c), and (d), indicate the points where 
the t w and t m distributions are shown in the corresponding figures (b), (c), and (d). The value of 
the load is 14, 1, and 0.1pN-nm-s-rad _1 , for (b), (c), and (d) respectively. Probability distributions 
P m for t m are shown by solid lines and P w for t w are shown by dashed lines. The waiting-times 
are exponentially distributed in all load range. The moving-times show peaked distributions. The 
peak at high load is partially hidden by a flat region for small t w , as a consequence of the rotor 
fluctuations. 
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FIG. 3: Dependence of the knee speed io n on r = k+/k- and on the cutoff 5 C for the symmetric 
potential with a parabolic bottom, (a) Motor torque as a function of the normalized angular speed 
for two different values of r: r = (black solid line), and r = 1.4 (red dashed line). Different from 
the V-shaped potential, there is a small plateau regime even for r = 0. (b) £ = u> n /uj 

max clS a 

function of r. (c) Torque-speed curves for different values of the cutoff 5 C . Black, red, green, blue, 
and violet lines correspond to values of 5 C of 1, 2, 3, 6, lOOJo respectively. Si and £2 correspond to 
the values of £ for S c = So and 5 C = W05o respectively, (d) £ as a function of the ratio 5 c /5o shown 
as the black solid line. The red dashed line corresponds to the asymptotic value £(<5 C — > 00). 
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FIG. 4: Torque-speed curves for an asymmetric potential. The torque-speed relationship is very 
similar to the symmetric potential case. The figure shows the motor torque r as a function of 
rotational speed for N= 1, 2, 8. For a given value of N the torque is almost constant up to the 
knee before it decreases roughly linearly. At low speed, the motor torque increases linearly with N. 
At high speed near zero load all curves collapse to the same maximum speed. These characteristics 
are apparent in the insert, where the torque per stator t/N versus tu are shown. 
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FIG. 5: Torque-speed curves for the symmetric potential with a smooth, parabolic bottom. The 
behavior is very similar to those observed in the symmetric and asymmetric V-shaped potential 
cases. The figure shows the motor torque r as a function of rotational speed for N= 1,2,..., 8. 
For a given value of N the torque is practically constant up to the knee, which is reached at 
almost the same speed value for different number of stators N, then it decreases linearly with u. 
The insert shows the "semi-parabolic" potential (solid red curve) and the symmetric V-shaped 
potential (dashed black curve) for comparison. 
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FIG. 6: Independence of the torque-speed characteristics on the load-rotor spring. Torques-speed 
curves (N = 8) with spring (dotted lines) and without spring (squared lines) show almost identical 
behavior. Instead, the concavity of the torque-speed curve depends on the stator jump rate ratio 
r = k+/k_: (a) r = 0, (b) r = 0.5, and (c) r = 1. 
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Quantity 


Value 


Comment 




0.002pN-nm-s-rad~ 1 


Estimated from [2] 


& 


« (0.002-50)pN-nm-s-rad- 1 






2vr/26 


From Ref. [6] 


^0 


(5/2 




T 


295.85K 


Room temperature 


TO 


505pN-nm 


Typical value 


T+ 


15KijT/0.95(5 


Typical value 


T~ 


15KijT/0.05(5 


Typical value 


fc + (symm. case) 


12000s- 1 


Fitting data 


/c_(symm. case) 


2fc+ 


From theory 


fc + (asymm. case) 


10000s" 1 


Fitting data 


L(asymm. case) 


20fc + 


From theory 


K 


400pN-nm-rad _1 


From Ref. [3] 


Kl 


4000pN-nm-rad~ 1 


From Ref. [3] 


(e - e L )o 


2vr/3 


From Ref. [3] 


A(0 - L ) 


2vr/7 


From Ref. [3] 


A£ 


^0 





TABLE I: Parameters used in the calculation. 
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